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ABSTRACT 

We  describe  a  framework  for  the  efficient  mod¬ 
eling  and  performance  evaluation  of  large  networks 
consisting  of  a  mixture  of  strategic  and  tactical  com¬ 
ponents.  The  method  emphasizes  hierarchical,  lay¬ 
ered  techniques  that  are  fed  parametric  models  at  the 
lower  level.  In  addition  to  the  algorithmic  structure, 
and  some  initial  algorithms  we  describe  an  object  ori¬ 
ented  software  architecture  that  is  under  development 
to  support  these  algorithmic  methods  in  a  distributed 
environment. 

1  INTRODUCTION 

Tactical  and  strategic  modeling  of  communication 
networks  is  a  daunting  challenge  given  the  complex 
nature  and  large  scope  of  today’s  and  tomorrow’s 
communication  systems.  The  approach  being  pursued 
within  ATIRP  Task  2.1  is  top-down  in  design,  but 
bottom-up  in  function-relying  on  low-level,  discrete 
event  simulation  of  tactical  communication  networks 
to  generate  statistical  parameters  for  higher-level,  hi¬ 
erarchical,  analytical  models  to  permit  fast  approxi¬ 
mate  performance  evaluation  of  strategic  communica¬ 
tion  networks. 

The  distinction  between  tactical  and  strategic  net¬ 
work  simulation  is  primarily  one  of  the  number  of 
nodes  to  be  modeled  in  the  simulation.  Tactical  net¬ 
works  are  classified  as  those  small  enough  to  be  simu¬ 
lated  in  reasonable  time  with  a  discrete  event  simula¬ 
tion  package.  All  other  networks  too  large  for  this  are 
designated  as  strategic,  and  must  be  evaluated  using 

‘This  work  was  supported  in  part  by  the  U.S.  Department 
of  the  Army,  Army  Research  Laboratory  under  Cooperative 
Agreement  DAAL01-96-2-0002,  Federated  Laboratory  ATIRP 
Consortium,  and  in  part  by  the  Center  for  Satellite  and  Hybrid 
Communication  Networks,  under  NASA  contract  NAGW-2777. 


approximation  methods. 

The  discrete  event  simulation  package  of  choice  for 
the  consortium  is  OPNET.  OPNET  will  be  used  to 
accurately  model  small  tactical  communication  net¬ 
works.  What  signifies  a  tactical  network  will  depend 
on  the  scenario.  One  may  consist  of  a  collection  of  mo¬ 
bile  SINCGARS-equipped  units,  or  perhaps  a  cluster 
of  mobile  users  surrounding  rapidly  deployable  cellu¬ 
lar  base  stations.  The  purpose  of  simulating  networks 
at  the  tactical  level  is  two-fold:  (i)  to  accurately  ver¬ 
ify  the  functionality  of  network  protocols  and  quantify 
their  operation  in  terms  of  desired  performance  met¬ 
rics  and  (ii)  to  provide  appropriate  statistical  param¬ 
eters  to  feed  into  the  lowest  level  of  a  hierarchical  an¬ 
alytic  model.  In  the  context  of  a  strategic  simulation, 
each  tactical  network  becomes  a  node  at  the  lowest 
level  of  the  mathematical  hierarchy.  These  nodes  are 
referred  to  as  “level  1”  nodes. 

The  strategic  network  performance  evaluation  tool 
is  currently  under  development  at  the  University  of 
Maryland.  The  tool  consists  of  an  algorithmic  struc¬ 
ture  and  a  software  architecture.  The  algorithm  is 
hierarchical  in  form,  deriving  its  structure  from  that 
of  the  network  which  follows  the  network’s  routing 
hierarchy.  The  mathematical  basis  for  the  model  is  a 
hierarchical  modeling  of  traffic  flows  based  on  Markov 
modulated  Poisson  and  Bernoulli  processes.  These 
models  have  been  known  in  practice  to  model  accu¬ 
rately  a  variety  of  multimedia  traffic  over  a  variety  of 
channels  and  communication  media.  At  each  level,  a 
node  represents  a  set  of  nodes  at  a  lower  level.  The 
boundary  of  a  node  is  characterized  by  a  set  of  sta¬ 
tistical  models.  The  statistics,  and  their  parametric 
representations  for  the  models  of  the  level  1  nodes  are 
obtained  from  the  OPNET  simulations  which  are  run 
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offline  . 

At  each  level,  topology  is  represented  by  logical 
nodes  and  logical  links.  At  the  lowest  (which  is  the 
highest  resolution)  level,  each  node  and  link  repre¬ 
sents  a  real  switching  system  and  a  physical  link.  At 
higher  layers,  each  node  may  represent  either  a  real 
system,  or  a  subnetwork,  or  a  group  of  switching  sys¬ 
tems.  At  higher  layers  links  may  correspond  to  either 
a  real  physical  link  or  a  virtual  link.  Nodes  are  col¬ 
lected  into  peer  groups  (clusters).  Peer  groups  are 
organized  into  a  hierarchy  where  one  or  more  peer 
groups  are  associated  with  a  parent  peer  group  (clus¬ 
ter)  and  so  on.  The  states  of  the  Hidden  Markov 
Models  that  we  have  developed  correspond  to  the  link 
states  as  they  are  used  in  hierarchical  state  dependent 
routing,  to  node  (switch)  states  as  they  are  used  in 
dynamic  adaptive  routing,  and  to  other  parameters 
of  the  multimedia  flow  including  Quality  of  Service 
(QoS)  parameters. 

The  model  defines  procedures  that  convey  model¬ 
ing,  simulation  and  performance  data  up  and  down 
the  hierarchy  of  models.  There  may  be  different  hi¬ 
erarchies  and  the  most  typical  one  is  based  on  the 
resolution  and  accuracy  (fine  or  coarse)  of  perfor¬ 
mance  and  traffic  data.  The  methodology  allows  user- 
definable  hierarchy  levels.  It  defines  how  each  level  of 
the  hierarchy  works  and  how  multiple  nodes  at  one 
level  can  be  summarized  in  the  higher  layer.  It  also 
defines  how  state  information  among  nodes  within  the 
same  level  and  between  different  levels  are  exchanged. 

The  models  in  this  approach  form  a  hierarchi¬ 
cal,  progressive  chain  (from  coarse  to  fine).  We 
will  then  develop  techniques  from  multi-level  opti¬ 
mization,  multi-grid  numerical  methods  for  progres¬ 
sive  computation  of  performance  metrics  or  trade-off 
curves.  Progressive  here  means  that  the  computations 
of  the  previous  layer  are  used  and  are  improved  upon 

fiiJm  A"n)  =  (t>{Jn+ lj  -^n+l)  +  An~*"  ^n+l) 

n  =  1,  •  •  ■ ,  N,  1  =  fine,  N  =  coarse.  Here  J  and  K 
are  two  performance  metrics.  The  subindex  indicates 
the  layer  or  accuracy  of  the  computation.  The  sec¬ 
ond  term  in  the  right  hand  side  indicates  the  correc¬ 
tion  term  that  utilizes  the  data  at  the  new  layer  and 
some  aggregate  variables  propagated  from  the  previ¬ 
ous  layer.  In  this  paper  we  describe  the  algorithms 
developed  todate  for  such  networks.  We  consider  the 
following  performance  metrics  only:  end-to-end  delay, 


throughout  blocking  probabilities,  cell  or  packet  loss 
rate. 

We  are  developing  fast  algorithms  to  compute 
trade-off  curves,  and  to  perform  interactive  trade-off 
analysis,  among  these  performance  metrics  by  link¬ 
ing  the  hierarchical  modeling  and  simulation  system 
to  multi-objective  optimization  packages.  In  this  way 
we  can  compute  performance  metric  sensitivities  with 
respect  to  network  parameters  using  the  simulation 
model;  a  particular  instance  of  this  idea  is  what  is 
called  simulation  differentiation. 

2  ALGORITHMS  BASED  ON  FIXED- 
POINT  METHODS 

The  first  analytical  method  that  we  have  inves¬ 
tigated,  is  that  of  fixed,  point  methods  for  estimat¬ 
ing  blocking  in  multihop,  multirate  integrated  ser¬ 
vice  networks  with  state  dependent  admission  con¬ 
trol  and  routing.  This  analytical  technique  is  at  the 
connection  level,  and  applies  to  synchronous  transfer 
mode  (circuit-switched)  services,  and  to  asynchronous 
transfer  mode  (ATM)  services.  These  methods  have 
been  shown  to  be  efficient  and  accurate  in  previous 
work  [3-5]. 

Recently,  fixed  point  methods  or  reduced  load  ap¬ 
proximations  for  network  analysis  have  attracted  sub¬ 
stantial  interest  [2-5].  They  provide  an  appealing  al¬ 
ternative  and  complement  to  discrete  event  simula¬ 
tions,  especially  in  network  design  and  in  performance 
evaluation  of  large  networks.  They  also  can  be  used 
for  worse-case  analysis  since  they  are  quite  accurate 
under  overload  conditions. 

In  our  work,  we  have  concentrated  on  integrated 
services  networks,  vs  voice  where  past  work  focused, 
and  on  hierarchical  solutions  and  formulations  of  the 
fixed  point  equations.  One  of  the  difficulties  of  moving 
these  effective  methods  into  multirate  methods  has 
been  the  fact  that  in  multirate  networks  the  compu¬ 
tational  complexity  of  these  algorithms  is  infeasible. 
Here  we  use  some  methods  first  introduced  in  [3]  to 
circumvent  this  computational  complexity  obstacle. 

The  network  model  that  we  assumed  is  as  follows. 

We  consider  general  networks  of  N  nodes,  indexed 
1,  2,  •  ■  ■ ,  N .  A  link  (i,j)  between  nodes  i  and  j  has 
capacity  C', j ,  counted  in  bandwidth  units  termed 
trunks.  Links  are  undirected;  (i,j)  and  (j.  i)  denote 
the  same  link.  Calls  offered  to  the  network  fall  into  S 
classes,  differentiated  by  bandwidth  and  possibly  by 
priority  via  trunk  reservation,  as  described  below.  A 
call  of  class  s  has  bandwidth  bs ,  meaning  that  if  the 
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call  is  admitted  to  the  network  then  it  is  allocated 
bs  trunks  on  each  link  of  a  path  from  its  source  to 
its  destination.  A  link’s  state  is  characterized  by  the 
number  of  calls  in  progress  ns .  of  each  call  class. 

A  call  is  admissible  on  a  given  path  if  it  is  admissi¬ 
ble  on  every  link  of  that  path.  Consider  a  call  of  class 
s  offered  to  node  pair  Link  admissibility  may 

be  decided  with  or  without  control.  In  the  absence 
of  control,  the  call  is  admissible  on  a  link  if  the  call’s 
bandwidth  bs  is  less  than  or  equal  to  the  idle  band¬ 
width  C  —  Admission  controls  considered  here 

are  of  trunk  reservation  type.  Specifically,  the  call  is 
admissible  on  link  (i,  j)  if  the  calls  bandwidth  bs  is 
less  than  or  equal  to  the  idle  bandwidth  C  —  T,tbtnt 
plus  r,  where  r  is  a  static  parameter,  which  may  de¬ 
pend  on  s,i,j ,  and  other  factors,  such  as  the  length 
of  the  route  associated  with  the  call.  Thus,  choosing 
parameter  r  to  be  small  allows  the  call  type  greater 
access  to  scarce  capacity.  If  admitted  on  a  path  in  the 
network  the  call  simultaneously  seizes  for  its  exclusive 
use  bs  trunks  on  every  link  of  the  path.  Similarly,  the 
call  simultaneously  releases  bs  trunks  on  every  link  of 
the  path  when  it  completes. 

We  consider  routing  policies  where  each  node  pair 
(i,  j)  is  assigned  an  ordered  list  Pi, 3  of  paths  between 

nodes  i  and  j:  P-1} ,  P;2} .  The  call  is  admitted  on 

the  first  path  in  the  list  where  it  is  admissible.  If  the 
call  is  not  admissible  on  any  path  then  it  is  blocked; 
that  is,  rejected  and  lost. 

Traffic  is  modeled  as  follows.  Calls  of  class  s  are 
assumed  to  arrive  to  node  pair  (i.  j)  according  to  a 
Poisson  process  with  rate  A ij(s).  If  admitted  to  the 
network,  the  call  holds  for  an  exponential  period  of 
time  with  mean  //,.,(>•).  Thus,  these  calls  constitute  a 
load  of  Pij(s)  =  Aij(s)/pij(s). 

The  key  performance  metrics  are  the  node  pair  by 
node  pair,  class  by  class  call  blocking  probabilities, 

Pr{blocking  a  call  of  class  s  arriving  to 

node  pair  (i,  j)}, 

or  (equivalently)  the  availability  probabilities,  defined 
as  corresponding  probabilities  of  admitting  the  call. 

The  fixed  point  method  that  we  use  is  based  on  the 
following  two  approximations  (assumptions): 

(Al)  (Link  Independence):  The  probability  that  a 
cell  is  admissible  on  a  given  path  in  the  network  is  the 
product  of  the  probabilities  that  the  call  is  admissible 
on  each  link  of  that  path. 


(A2)  (Poisson  Rates):  Arrivals  of  calls  of  a  given 
class  to  a  given  link  are  described  by  a  Poisson  pro¬ 
cess.  The  rate  of  the  process  may  depend  on  the  state 
of  the  link. 

It  is  important  to  note  that  although  arrival  rates 
to  node  pairs  are  given  and  are  assumed  to  be  Poisson, 
admission  control  and  routing  depends  on  link-states 
and  thus  perturbs  the  Poisson  character  of  the  link 
arrivals.  These  two  assumptions,  lead  to  a  system  of 
equations,  whose  unknowns  are,  for  each  link  ( i,j ) 
and  call  class  s: 

(i)  arrival  rates  i gj(s)  of  calls  of  class  s  from  any 
node  pair  that  includes  link  (i,  j)  on  some  route, 
given  that  link  (i,  j)  is  in  a  state  that  admit  calls 
of  class  s,  and 

(ii)  probabilities  a{j(s)  that  link  (i.  j)  is  in  a  state 
that  admits  calls  of  class  s. 

The  fixed  point  formulation,  essentially  consists  of 
two  mappings.  One  such  mapping  is  for  every  node 
pair  (i.  j)  from  probabilities  akj(s)  (with  k  and  l  de¬ 
pendent  on  the  routing  structure)  to  rates  i 'i  f(s).  The 
other  mapping  is  for  each  link  ( i,j )  from  rates  Vij{s) 
to  probabilities  oy(s). 

To  estimate  equilibrium  network  performance,  we 
solve  for  the  fixed  point  of  the  composition  of  these 
two  mappings,  simultaneously  fixing  the  Pi,j{s)  and 
aij(s).  The  algorithms  as  in  [3]  are  inspired  by  the 
close  mathematical  similarity  of  these  equations  to 
certain  problems  in  network  reliability. 

The  mapping  from  admissibility  probabilities  to  ar¬ 
rival  rate  is  now  described. 

Let  v\™\k,l,s)  =  rate  that  ( k,l )  admits  class  .s 
calls  on  its  rnlh  route,  if  that  route  contains  link  (i,  j); 
otherwise  (k,l,  s)  =  0.  Thus,  link  (i,j)  receives 
calls  of  class  s  at  aggregate  rate 

Vi,j{s)  =  T,mT,k^j){k,l,s). 

We  now  proceed  to  compute  v^™\k,  l,  s),  assuming 
link  (*,  j)  belongs  to  the  rnlh  route,  for  node  pair  ( k ,  l). 

For  m  =  1 

^  s)  =  h,lis)  II  au,v{s) 

where  A k,i(s)  denotes  the  original  rate  of  arrivals 
for  calls  of  class  s  for  node  pair  (k,l)  and  aUjV(s)  is 
the  admissibility  probability  for  class  5  on  link  {u,v). 
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For  m  >  1, 

^f(k,l,s)  =  A  kj(s)-R%](k,l,sy 

n (u,v)eP^\(u,v)^j)  au ^ 

R^j\k,l,s)  =  Pr {none  of  the  paths 

P^1}  ■  ■  ■  p}™'  i')  are  admissible  I 
path  P-™  is  admissible} . 

Combining  the  formulas  above  completes  the  map¬ 
ping  from  link  admissibility  probabilities  to  link  ar¬ 
rival  rates. 

We  next  describe  the  mapping  from  arrival  rates 
to  admissibility  probabilities  when  there  is  admission 
control  in  the  network.  It  is  well  known  that  admis¬ 
sion  control  destroys  the  product-form  of  the  link  oc¬ 
cupancy  probabilities  P(n ),  which  in  turn  destroys  the 
efficient  and  exact  computation  of  those  probabilities 
by  a  formula  such  as: 

"(*)  =  1  -  >:u  a-  b,  •  i  /’(») 

where  C  is  the  link  capacity,  S  the  number  of  classes 
and  bs  the  bandwidth  of  each  class  s  call.  Instead, 
we  now  have  to  compute  the  equilibrium  distribution 
of  a  Markov  chain  with  a  large  state  space;  a  pro¬ 
hibitive  computational  task.  Instead  we  employ  an 
approximation  first  introduced  in  [6,7].  This  method 
transforms  the  problem  into  one-dimensional  one  es¬ 
sentially  by: 

(a)  assuming  that  while  ns ,  the  number  of  calls  in 
progress  of  a  class  s  varies,  the  proportion  of  such 
calls  in  progress  remains  fixed  (or  varies  slowly), 
and 

(b)  obtaining  that  proportion  by  equating  the  admis¬ 
sion  rate  of  calls  of  class  s  to  the  corresponding 
departure  rate. 

The  following  method  appears  to  be  very  effective. 
Letting  a(s)  denote  the  average  number  of  calls  of 
type  5  in  progress, 

a(s)  =  a(s)\(s)/ p,(s). 

since  these  calls  enter  into  service  at  rate  a(s)A(s)  and 
depart  at  rate  a(s)/u(s).  Treat  the  a(s)  as  given,  and 
consider  the  one-dimensional  Markov  chain,  which, 
for  any  state  n  and  call  classes  s,  jumps  to: 


(c)  state  n  +  bs  with  rate  A (s)I(C-n  >  r  +  6s),  where 
r  is  the  trunk  reservation  parameter,  and  to 

(d)  state  n  —  bs  with  rate  p(s)n^^I(n  >  bs).  In 
the  true  system,  this  jump  down  depends  on  the 
number  of  calls  in  progress  nt  for  each  class  t, 
rather  than  the  aggregate  state  n  =  nt ,  so  this 
is  an  approximation,  that  follows  from  estimating 
ns  by  na{s)/T,ta(t). 

Then  the  probability  of  admitting  a  call  of  class  5 
is  approximately  a(s)  =  1  —  H^=c_bs+1P(n).  Tak¬ 
ing  P(n )  to  be  the  equilibrium  distribution  for  this 
one-dimensional  chain,  we  arrive  at  a  small  system  of 
fixed  point  equations.  In  [6,7]  it  is  shown  that  this 
approximation  becomes  asymptotically  exact  in  over¬ 
load  conditions,  for  large  C  and  bs  «  C.  More  re¬ 
fined  approximations  are  also  discussed  in  [6,7]. 

In  our  work  todate  we  have  developed  hierarchical 
versions  of  these  algorithms  that  are  based  on  hier¬ 
archical  routing  in  the  network  being  analyzed.  We 
are  currently  working  on  versions  of  these  hierarchical 
algorithms  that  are  progressive. 

3  SOFTWARE  ARCHITECTURE 

The  software  implementation  of  the  performance 
evaluation  tool  is  written  in  Java.  Consortium  users 
will  access  the  tool  via  any  World  Wide  Web  browser 
that  supports  Java  from  any  hardware  platform.  The 
software  architecture  is  fundamentally  client/server- 
based.  Users  access  the  simulation  kernel-hosted 
on  a  server  at  the  University-via  a  browser-based, 
graphical  user  interface  (GUI)  which  forms  the  client. 
Thus,  while  simulations  are  constructed  within  the 
client,  they  are  executed  within  a  high-speed,  multi¬ 
processor  server. 

The  software  architecture  is  CORBA-compliant, 
with  all  Java  objects  being  specified  first  in  CORBA’s 
Interface  Definition  Language  (IDL).  The  Object  Re¬ 
quest  Broker  of  choice  is  Iona’s  OrbixWeb.  Simulation 
objects  may  be  made  persistent  via  an  object  adapter 
to  ObjectStore. 

Nodes  and  channels  form  the  basis  of  the  network 
model’s  object  structure.  Nodes  form  an  object  hier¬ 
archy  of  their  own,  and  have  a  different  representation 
on  both  the  client  and  server.  Server  node  objects 
contain  all  data  that  must  be  persistent  from  run  to 
run,  and  client  node  objects  are  decorated  versions 
of  the  server  objects  containing  additional  GUI  func¬ 
tionality.  Channel  objects  also  have  differing  repre¬ 
sentations  on  both  client  and  server.  Client  channel 
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Figure  1:  The  GUI  to  the  hierarchical  performance  evaluation  tool. 


objects  primarily  capture  the  spatial,  topological  and 
qualitative  characteristics  of  the  channel  whereas,  on 
the  server,  channel  objects  are  translated  into  math¬ 
ematical  quantities  imbedded  in  the  model. 

Figure  1,  depicts  the  current  GUI  to  the  hierarchi¬ 
cal  performance  evaluation  tool. 

The  near-term  goals  of  our  work  are  to  develop  a 
Java-based  GUI  prototype,  and  to  have  a  preliminary 
alpha  version  incorporating  a  single  “level  1”  node 
type  for  mathematical  analysis  by  August  1997.  Thus, 
the  initial  version  will  be  capable  of  approximate  per¬ 
formance  analysis  over  a  large  homogeneous  network. 
Introducing  heterogeneity  into  the  system  will  be  the 
subject  of  future  work. 
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